Generating multi-scaling networks with different types of nodes 
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A variety of scale- free networks have been created since the pioneer work by A.-L. Barabasi and 
R. Albert. All this networks are homogeneous since they are composed of the same kind of nodes. 
In the realistic world, however, one element (node or vertex) in the network may play different roles 
and hence has different functions. In this Letter, we develop a new kind of network to account for 
this property. In our model, each type of nodes may exhibit a scaling law in the degree distribution 
and the scaling exponents are adjustable. As a consequence, the whole network lacks of such scaling 
characteristics, which indicates that many previous statistical results based on empirical data that 
claimed to be scale-free networks may need to be reexamined. This model poses an alternative 
way of the network division other than the module method. Besides, one can expect that this new 
network will exhibit some interesting properties concerning the dynamical processes on it. 

PACS numbers: 89.75.Da, 89.75Fb, 89.75Hc. 



In recent years, an increasing interest on complex net- 
works has risen in physical society because it shows 
as one of the most promising tools to describe various 
social, biological, and technical systems, such as the 
Internet H ll H, the World Wide Web (WWW) 4] or 
collaboration networks 5, 6J. In the models, vertices or 
nodes may represent people, proteins, species, routers, 
or html documents while the links between the nodes 
correspond to acquaintance, physical interactions, preda- 
tion relationships, cable connections, or hyper- links, re- 
spectively. In real- world networks, two outstanding phe- 
nomenon are mostly presented. One is the small-world 
phenomenon 0,13) which refers to that one can promptly 
reach remote parts of the network through very few hops. 
The other is the scale- free phenomenon 9^ ilOL ilL] , which 
refers to the degree of the nodes (the number of links 
connected to it) obeying a power-law distribution as 
P{k) - k-'i. 

To explain the scale-free property of a network, A.- 
L. Barabasi and R. Albert (BA)'9'| designed a model by 
introducing the concepts of growing network and of pref- 
erential attachment. The network are growing at a con- 
stant rate and new nodes are attached to the older ones 
with a probability which is proportional to the degrees 
of the preexisting nodes. 
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The network generated in this way has a fixed scaling 
exponent 7 = 3. Variant models that assign a fitness to 
every newly-added node, which accounts for the impor- 
tance of the node for attracting new links, can gen erate 
networks with adjustable scaHng exponents [l^ Tl^ . 
Besides, many other network models and methods are 
proposed to probe the various properties in real-world 
networks El [H [H [11 El . 

Nevertheless, all the previous works are focused on net- 
works in which all the nodes are identical and play a 
unique function in the homogeneous (or inhomogeneous) 



network. As one explicitly examines into a real-world 
network, however, it can be found that a node is gener- 
ally not play only one role in the network. It may play 
multi-roles or couple with nodes in other systems. As 
show in Fig.l, nodes of type A or type B respectively 
form two subsystems with links between some of them. 
From another point of view, if we regard the whole object 
as two different networks which are respectively consisted 
of type A or type B nodes, then a node in one network 
may couple with a node in the other network even though 
this two networks may have completely different topology 
or function. These nodes may behavior quite differently 
as they present in different environments. For example, 
a researcher may work in two or more fields and coop- 
erate with different authors. He writes papers in con- 
densed matter physics with collaborators who are always 
engaged in this field while he may also contribute to the 
high energy physics by participating another group. He 
may also have an independent friendship network |20j. 

Many works have been carried out to deal with the 
community or module structure of the network, where 
nodes is liable to interconnect within one community 1 1 (t 
nil m . In this Letter, we introduce another way to 
catalogue the network nodes according to their functions 
as relative independent subsystems. We will construct a 
heterogeneous network consisted of two types of nodes, 
type A and type B. As in a biochemical reaction pro- 
cesses, the reactants and their relationship form a com- 
plex network. One type of reactant preferentially com- 
bines with some other specific reactants. The attractance 
or binding energy ai between the same type of nodes (A- 
A or B-B) differs from the attractance a2 between the 
different types of nodes (A-B). When a new node (type 
A or type B) joins the growing network, two rules govern 
its attachment to the preexisting nodes: i) The Barabasi 
and Albert preferential attachment, and ii) The selection 
of partners according to its genus. The whole network is 
such composed of two types of nodes or vertices, i.e., it is 
a heterogeneous network with two connected subgraphs 
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of Ga and Gb- Figure 1 is a sample network generated 
by this algorithm for 30 nodes with 02/ ai = 0.1 and the 
proportion of type A nodes pa = 0.6. 

It is found that our model network exhibits a multi- 
scaling structure in which the degree distribution for each 
subgraph Ga or Gb shows a scaling law. As a conse- 
quence, the total degree distribution lacks of such scaling 
characteristics except for some special cases. The result 
is instructive for many recent statistics on empirical data 
that claimed to have the scaling law distribution. They 
are generally not reliable if the nodes play multi-roles 
in the generating process or the network is composed of 



several subsystems and need to be carefully reexamined. 
Many of them may show as the pseudo-scaling laws. 

We assume that the newly-added node presents as type 
A with probability pA and as type B with probability 
= 1 — PA, respectively. Each node has m feet to 
be connected to the existing network. The vertices of 
type A and type B respectively form two vertex sets of 

Va = {vi,V2, - ■ ■ ,vna} and Wb = {wi,W2, ■ ■ ■ ,wnb}- 
Here Na and Nb are the total numbers of node A and 
node B, respectively. From the above rules, we obtain 
the continuously growing equations 



dkAji) _ mpAaikAji) ^ mpBa2kA{i) 

dkBji) _ mpAa2kB (») ^ mpBaikB (i) 



where fcyi(i) and kB{i) are the degree of the i-th node 
of either type A or type B, respectively. The first line 
of the above equations describes the growth rate of type 
A nodes in the network. The first term on the right- 
handed side represents a newly-added type A node, which 
is generated with probability pA, is attached to a type A 
node in the preexisting network while the second term 
is a newly-added type B node to be attached also to a 
type A node in the network. Analogously, the second line 
describes the growth rate of type B nodes in the network. 

To solve the above equations, we consider the ther- 
modynamical approximation. Suppose the system has 
multi-scaling law, namely, the dynamic exponents de- 
pend on the attractance ai and a2, 

kA{t,ta,t) = m{-^f^ 
to 

kB{t,to,i) = m{^f^, (3) 



where Iq is the time at which the node i was born. The 
dynamic exponents Pa and (3b are bounded, i.e. < 
PaiPb < 1 because a node always increases the number 
of links in time {f^A, Pb > 0) and kA,kB cannot increase 
faster than t (Pa,(3b < !)• We calculate the sum over 
j in Eqs.(|2Jl by writing them in the integral forms. By 
noting that the nodes respectively belong to sets Va and 
Wb with probability pA and p b , one has 

j:^kAij) = mJ^PAi-r^ - 

Here we have dropped the term for it becomes less 
important as < — > 00. 

After some calculations by introducing a variable z — 
(1 — — Pb), we obtain a third order equation as 



I 

z3 + [^(^-,^^-l) + (^^-2)]z^ + [^(-^-^^ + l) + 4(2-^)].-4=0, (5) 
Pb 012 ai ai pb 012 ai p% ai 



and 



PAOLl+PBa2Z PBaiZ+pAa2 

Therefore, it demonstrates that in the thermodynamical 
limit, the degree distributions of the constituent nodes 



respectively obey the power law. The power exponents of 
each type of nodes are generally different. Obviously, the 
total degree distribution that comprises both type A and 
type B nodes will not exhibit the scaling characteristics. 

Figure 2 plots the dependence of the exponents ja and 
7b (7 = 1 + on the relative attractive strength 
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between the nodes aijax. It is seen that the expo- 
nents can be either larger than 3 or in the regime of 
2 < 7 < 3, depending on the parameter value of ailoL\. 
The upper panel shows that if 7^1 > 3, then 7^ < 3 
whereas if < 3, then 75 > 3. For smaller a2/ai, 
the newly-added node will preferentially attaches to the 
same type of preexisting nodes. In the case of > 0.5, 
type A nodes with larger degrees are continuing to at- 
tract more links and thus the exponent is relative small. 
As aila\ — > 0, the scaling exponents for both genuses 
of nodes tend to 3. In this case, the coupling between 
the two types of nodes are weak and the resultant net- 
work nearly divides into two isolated subsystems. At the 
point of aij a\ = 1, which corresponds to the case that A 
and B nodes are completely identical, our model degen- 
erates into the ordinary BA model with a fixed exponent 
7 = 3. As a special case of j)a = Vb = 0.5, the exponents 
= 7b = 3, regardless of the detailed values of a^joti- 
For larger 02 /ai (a2/ai > 1, see the lower panel of 
Fig. 2) where one type of node will preferentially attach 
to the different type of nodes (A to B or B to A), the 
network seems to be interweaved alternatively by node 
A and node B. Hence for ipA > 0.5, it becomes hard 
for the type A nodes that already have more links to 
get even more links, i.e., the degree growth is damped. 
Most nodes of this type will remain few links and so the 
scaling exponent becomes larger. On the other hand, 
for pb < 0.5, those nodes of type B that already have 
more links will continue to attract more links from type 
A nodes and so the scaling exponent decreases as ctijcti 
increases. 

To check the above theoretical results, we simulate the 
growth process on the computer. We start with Nq — \ 
interconnected nodes consisted of 2 type A and 2 type 
B nodes, respectively. Figure 3 displays the simulated 
degree distribution P(fc) for pA = 0.9 and a2/ai = 0.5. 
A total of 200,000 nodes are involved. The solid lines are 
respectively the theoretical results for type A (7^ = 2.95) 
and type B (75 = 4.06) nodes. It is seen that the simu- 
lating data coincide with the theoretical result quite well. 
Here we point out an important fact. Just from the sim- 
ulating data (see the inset of Fig. 3), one may misinter- 
pret that the total distribution indiscriminating the node 
genuses also show a power law relation. However, it is 
incorrect from the theoretical considerations. This fact 
precautions us that it should be wary when dealing with 
empirical data that include several types of nodes to draw 
out a scaling law. 

Figure 4 is the same as in Fig. 3 for 02/011 = 3.0. The 
two solid lines correspond to the theoretical results with 
7A — 3.96 and 7s ~ 2.18, respectively. It shows that for 
type B nodes, the data deviate form the theoretical lines. 
This is because we have dropped in Eqs.Q) a term which 
behaviors as . When /? ^ 1 (or 7 2), as in the 
present case, this term becomes increasingly important 
for limited number of simulating nodes and hence the 



approximation becomes poor. With the increase of node 
number, the coincidence should improve correspondingly. 

In summary, we have developed a bipartite network in 
which the nodes are divided into two genuses according 
to the interactions between them. There is an important 
difference between our model and the fitness model. In 
the fitness model, each node has fixed fitness while in our 
model, the interaction between a node with others is de- 
pendent on the types of its partners. Just as in a library, 
one classifies the books into catalogues and the books are 
connected by the cross-index table. One can reach a spe- 
cific book through different ways by following the classifi- 
cation method. Two major conclusions can be reached in 
our model: i) The network is catalogued by the genuses 
or functions of the nodes while most previous works di- 
vide the network by communities or modules. Our model 
is in fact heterogenous, ii) There are multi-scaling char- 
acteristics for each type of nodes, which implies that the 
total degree lacks of a power law distribution and many 
previous empirical statistics may need to be reexamined. 
Our model may provide a prototype to discuss couplings 
between two or even more subsystems which have scaling 
properties in degree distributions. It is also interesting 
to explore various dynamical processes such as searching 
processes 22J on this network, we expect that new algo- 
rithms based on our model or its possible variants will 
largely promote the searching efficiency on this kind of 
networks. 
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Figure Captions 
Figurel (Color online) A heterogeneous network gen- 
erated from the rules described in the context. The dot- 



ted lines represent interactions between different types of 
nodes. 

Figure 2 (Color online) Dependence of the scaling ex- 
ponents of the degree distributions for either type A (7^) 
or type B (75) nodes on parameter ratio a^joii- Upper 
panel: = 0.9. Lower panel: 7^ versus Oi^jotx for 
■pA = 0.05, 0.20, 0.50, 0.70, 0.90. 

Figure 3 (Color online) Comparison of simulating de- 
gree distributions for type A and type B nodes with the- 
oretical results. a2/ai — 0.5. The total number of nodes 
are 200,000 and pa = 0.9. Inset: The degree distribution 
indiscriminating the node genus. 

Figure 4 (Color online) Same as in figure 3 for 0:2/01 = 
3.0. 



